Analytical Current Formula for Strongly Coupled Quantum Dot with Arbitrary 

Electron-Electron Interaction 
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An analytical expression for the current through a single level quantum dot for arbitrary strength 
of the on-site electron-electron interaction is derived beyond standard mean-field theory. By describ- 
ing the localised states in terms of many-body operators, the employed diagrammatic technique for 
strong coupling enables inclusion of electron correlation effects into the description of the local 
dynamics, which provides transport properties that are consistent with recent experimental data. 



PACS numbers: 72.10.Bg, 73.63.Kv, 73.23.-b, 71.27. +a 



c3 



-a 
c 

o 
o 



> 

On 
OO 
(N 
(N 
O 

in 
o 
-i— > 



o 
o 



X 



Electron transport through few electron islands is 
strongly influenced by electron correlations between its 
atomic like states. A variety of correlation effects, in- 
cluding Kondo effect Q and resonant current peaks Q, 
have been observed. For low conductance of the tunnel 
junctions one can employ "orthodox theory" [3j, which 
treats the tunnelling in lowest order perturbation the- 
ory (golden rule). In the strong coupling regime, how- 
ever, this theory collapse since then the transport is not 
dominated by sequential tunnelling. While effects of co- 
herent tunnelling have been extensively studied Q, the 
question of a general description which includes correla- 
tion effects and is valid both for strong and weak cou- 
pling regimes remains open. In this Letter, an analytical 
formula for the current through a single level quantum 
dot (QD) is derived for arbitrary on-site Coulomb repul- 
sion, beyond standard mean-field theory (e.g. beyond 
self-consistent Hartree-Fock, or Hubbard I, approxima- 
tion (HIA) HIE0)- 

This formula include local correla- 
tion effects which provide the well known renormalisation 
of the localised level 0,Eli|- 

The present derivation concerns single level QDs, here 
described by the Anderson model e -g- 

« = E e ^ c U+ E E P hP 

kaeL.R P=0,<r,2 

+ E k<i(^+^ M )+M (i) 

kaeL,B. 

where c} ka (k) creates (annihilates) an electron in the left 
(L) or right (i?) lead at the energy and spin a =\, |, 
whereas X™ = \p)(a\, \q) = |0), \a), |2) = | U), is a 
Hubbard operator |5(, where h p = X pp is the population 
operator for the state at the energy E p . Here, the local 
state energies are Eq = 0, E a = e CT , and Ei = e-\ + 
ei + U, where e a is the localised single electron energy, 
whereas U is the on-site Coulomb repulsion. The last 
term in Eq. JQ) models the tunnelling, with rate Vka, 
between the leads and QD. The factor rj a = (a\d a \2) 
(rj^i = ±1) accounts for the selection rule between the 
singly and doubly occupied states ((Ole^ltr) = I), and a 



is the opposite spin of a. 

Using the description of the Anderson model in terms 
of Hubbard operators is motivated because this enables 
the use of a diagrammatic technique which i) provides an 
expansion of the QD GF in the strong coupling regime 
[1 1| . ii) yields the logarithmic renormalisation of the lo- 
calised level E2> an d Hi) can be employed for a direct 
transformation of the on-site properties into correspond- 
ing non-equilibrium quantities |13t 

The stationary current through the system can be cal- 
culated from the general expression Eq. (19) in Ref. 
with the matrix T L / R (given later) describing the cou- 
pling between the left /right lead and the QD, and the 
retarded/advanced/lesser matrix GF G r / a / < for the lo- 
cal QD properties. 

The components of the GF are defined by 0, flil H^ . 
El G o5 (f,t') = (-i)(TSX a (t)X l (t'))/(TS), where a, b 
denote Fermi like (single electron) transitions whereas 
a, b denote the conjugate transitions. The action S = 
exp(-tJl B -^H'(t)dt), 1/(3 = k B T, [H El EEl 03 
contains time dependent disturbance source fields U^(t), 
used for the generation of a diagrammatic expansion 
of the GF in terms of functional derivatives. Here, 
H'(t) = U Q (t)h° + J2 aa , XJ aiy ,{t)X™' + U 2 (t)h 2 . Phys- 
ical quantities are accquired from the GF as U^(t) — * 0. 

In the present system, there are four possible single 
electron transitions, e.g. X 0a and X a2 , which yields a 
4x4 matrix equation for the GF. Following the derivation 
in Refs. |l2|,[l_3|, one finds that the general structure of 
the equation of motion for the matrix GF is given by 



i^I - A - U(t)) G(t, t') = 6(t - t')P(t) 



-[P(t+) + R(i + )] 



-i/3 



V(t,t")G(t",t')dt", (2) 



where I is the identity, A = diag{A^ , Aj , A° x , A^ T } 
(diagonal matrix) contains the bare transition energies 
A° = E a - E and Af^. = E 2 - E s . The interaction 
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between the leads and the QD is given by 

w(t t>\-( V( *' t ' ) v '(*'*') 
v^,tj- ^v'(t,f) V(t,f) 

where V = diag{Pi-, Vj_} and V = cr 2 V (a z is the z- 
component of the Pauli spin vector), with V a (t,t') — 
Y,keL,R \ v ka\ 2 gka{t,t'). Here, gka(t,t') is the GF for the 
conduction electrons in the leads, satisfying the equation 
(id/dt-e ka )g ka (t,t , )=8(t-t'). 

In Eq. J21 the end-factor P(t), which arise due to 
the non-commutative algebra of Hubbard operators, car- 
ries the spectral weights of the components of the ma- 
trix GF. In the considered system P = diag{Pi,P 2 }, 
where P„, n = 1,2, themselves are 2x2 matrices 
given by P 1 = {PcwoW = {(T{X° CT , X*' }^}^, 
and P 2 = {P S2 2 S >}™> = {<T{X ff2 ,X 2ff '}) a } CTCT , (sub- 
script U signifying the source fields dependence of 
the averages). In the present case PoaaQ, Ps22a = 
0. The non-commutativity of the Hubbard opera- 
tors also gives rise to the operator R(i) 0, 0, Il3| . 
here given by R = diag{Ri,R 2 }, where Ri(£) = 
{Ro**'o(t)}a*' = {8 aa >8/8U Q {t) + 8/8U a . a (t)} aa , and 
Ra(t) = {R S 22a>{i)}„<r> = {8 SS . 8 / 8U 2 {t)+8 / 8U SS , {t)} aa , 
that is, functional differentiation operators which operate 
on the GF G with respect to the auxiliary source helds 
U c (t). The matrix U = diag{Ui,U 2 }, where Ui(t) = 
{U aa ,{t) ~ 8 aa ,U (t)}™' and U 2 (t) = {& 33 >U 2 (t) - 

The structure of the equation of motion, Eq. © , sug- 
gests that G(M') = D(i, i')P(t'), where D(t,i') is the 
locator \Ti\ carrying the local on-site properties of the 
GF, e.g. the position and width of the localised level. 
Here, fluctuations of the spectral weights are neglected, 
hence, P is treated as a constant with respect to the 
source fields U^(t). Therefore, varying the GF with re- 
spect to the source fields generates 8G = (<5D)P, and 
since the locator satisfies the matrix property DD 1 = 
/ = D *D one finds that 81) = -D((5D" 1 )D, hence 
SG = — D(JD -1 )G. Thus, it is necessary to study the 
locator and its inverse. 

In general, the equation of motion for the locator ap- 
pear very similar to Eq. that is, replacing G by D 
and removing P in the first term on the right hand side of 
Eq. (J2J. In addition, the structure of the resulting equa- 
tion of motion for D suggests that the inverted locator 
D 1 can be written 

D -1 ^,* 7 ) = d" 17 *,*') - S(t,t'), 

where the self operator S(t,t') [Tl| is identified by 



S(t,t') = {[P(t+) + K(t+)} 



to-ip 



V(t,t 1 )B(t 1 ,t 2 )} 
xT>^ 1 (t2,t')dt 2 dt 1 , 



present derivation, the dressed inverted locator in SG = 
— D((5D _1 )G is replaced by its corresponding bare quan- 
tity, e.g. 8G = — D(<5d _1 )G, meaning that the diagram- 
matic expansion terminate after the first functional dif- 
ferentiation, the so-called loop correction [ill Il2t Il3| . 
This will then generate the renormalisation of the lo- 
calised level, as is shown below. However, continuing 
the functional differentiation to higher orders generate 
higher order diagrams that account for additional many- 
body correlation effects [TlL Hjjj . for instance contribu- 
tions from the Kondo effect. 

Due to the non-commutativity of the matrices in Eq. 
it is now preferable to study the scalar equations. 
The interaction matrix V goes unaffected by the func- 
tional derivatives, since here the source fields are only 
applied to the localised QD states. Noticing that, for in- 
stance, R ^a(t+)d- 1 {t,t') = -8{t - i')iWo(i + )U(i), 
one finds that the only non-vanishing contributions are 

+D a2 „(t",t+)G^ l (t+,t')]. 

The Fourier transformed components Goo-a-ch G s2 a0 of 
Eq. (J2J, hence, reduces to 

(ilV - A ct0 - PQaaoV a )Goaao(i^) = PoaaO 

+ [V<yPoaaoV a + "t] s 8 A 2ff ] G s2 a0 (iCt>) , 

n ,. \ r\aPs22aVa + T]<j8A<j 

&ff 2 <r0 = ^-G 0aC ro{lUJ), 

VjJ - A 2a - Pa22a V a 

(as U^(t) — ► 0) and analogously for the components 
Ga22s, G s2 ao, with A a = + 8A a and 



"—5 E £ 

a—L.R k£a 



f a {e k a) - /H 

Ska - OJ 



x {-Zbn[Do*M + V*D r o*n(u)]}d»o, (3a) 



T O U r- ~. J 



f a {e k s) - /H 

Ska — W 



x {-2lTn[Dl s2a {u)+^Dl 22a {uj)]}dw, (3b) 

where D r r{u)) is the retarded form of the corresponding 
locator, and /l/.r(w) = fi^ — ^L/R) is the Fermi function 
at the chemical potential Hl/r for the left/right lead. 
The renormalisation due to the loop correction and its 
effects on the transport properties is discussed elsewhere 

The formal structure of the equation for the QD GF, 
e.g. Dyson equation, gives by a direct application of 
the Langreth rules for analytical continuation ^3 the 
retarded/ advanced and lesser counterparts of the GF as 



whereas the bare locator d(t, t') satisfies the equation 
[id/dtl - A - V(t)]d(t,t') = 8{t - t')I. For the 



G r / a \uj) 
G<\uj) = G r (uj)V<{Lu)G a (uj 
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Here g is the GF for the QD in the atomic limit, 
whereas the components of V < (oj) are given by Vf- (ui) = 

i[/i(w)rj + f R (u>)rf}; r CT = r£ + r R and v L J R = 

— 2Im^ fcgi ^ i? l^fecr 1 2 fffccr ( C<J ) ■ The retarded/advanced and 
lesser GFs can now be inserted into the general expres- 
sion for the current, Eq. (19) in Ref. (lif . where the trace 
is taken over the 4x4 matrices. It is then straightforward 
to show that 

tr [f L {oj)V L - f R (u)T R ][G r (u) - G»] = 

= -»X>4")r£ - f R (u)r R ]r a \gr(u)\\ (5a) 

tr [T L -T r }G<(uj) = 



T L a T R \[f L (u)T L a + /«Mr«]|^M| 2 ,(5b) 



with g r a = g\ 



OaaO 



vAGo^a + G l 2rT0 ] + G s22tf , giving 



Q r a (u) = (uj-A° a0 -SA a0 -5A 2& -UP 0a ) 

x[(w - A° - V7)(w - A° 5 - <5A ct0 - 8A 2S ) 
-U{P 92 VZ-8A aQ )]-\ (6) 

where P 0(T = Poaao and P ff2 = P&22&, and F CT r = V^(w). 
Hence, the current becomes 



r^|^M| 2 [/LM-/KM]^. (7) 



It should be emphasised that the derivation of Eq. J7J is 
valid for arbitrary couplings T L / R , not only T L = \T R . 

Here, the end-factors are given by Po<r = 
Im/[G<» - G>H - E CT 'G >,(u;)]da;/(27r) and 
P ff2 - Im/[G < ( w ) - G> ( W ) + E < G<, a H]dw/(27r) 1 
which leads to the boundary condition Po CT + P 5 2 = 
— tr Im J" G r (w)ciw/7r = 1 [l8| by the requirement that 
the integrated total local density of states is unity 
[ill Hence, Eqs. © and QJ, along with of 

the end-factors, lead to a set of equations subject to 
self-consistent calculations, for each value of the bias 
voltage. This, then provide an accurate field dependence 
of the transport properties of the system. Replacing 
T L - T R by T and f L T L - f R T R by f L T L + f R T R in Eq. 

one easily shows charge and current conservation of 
the presented result, e.g. d/dt(P 0a + P s2 ) = 0. 

Now, since the occupation number (n a ) = (h a + h 2 ) = 
Po-2, it follows that Pq„ = 1 — (riff). Hence, removing the 
renormalisation energies in Eq. ©, reduces the GF to 
the standard mean-field result, e.g. the HIA 0. 

The atomic limit of Eq. © is trivial since this limit 
is the starting point of the expansion of the QD GF. 
However, it is worth to notice that 5A a Q,SA 2s — > as 
Vka — > 0, c.f. Eq. 10. Further, in the limit of no on-site 
Coulomb interaction, e.g. U — > 0, the renormalisation 
energies SA a0 + SA 2& — >• 0, giving Q^{uj) = (u) — eq- 
VJ) _1 , i.e. the result from the exactly solvable Fano- 
Anderson model 0, H(| . Inserting this expression into 



the current, Eq. J7J, yields the well-known result for 
the single resonant level 0, 0|. Finally, in the limit 
of strong on-site Coulomb repulsion, e.g. U — > 00, Eq. 
© becomes Q^(uj) = P 0ct /(^ - A ct0 - PqctVJ), where 
A aQ = A° + 6A aQ and (at T = K) 



(n a ) 



Poa 

T7 



E r " 



a=L, R 



rrO 



— arctan ^ 

Tr P 0CT r ff /2 



since 5A rr o J 5A 2& are finite for all U. This GF was derived 
in Ref. 01 and is consistent with the result in Ref. 
By removing the loop correction SA a0 , this GF reduces 
to the result in Ref. 0. Also in this case the current 
becomes particularly simple, e.g. 



J 



E 



rirfl p2 



- A„ a - P aa VX\ 2 



[/z-M - f R (uj)]duj. 



The difference between this expression and the one for 
the single resonant level is that the end-factor P$ a makes 
sure that the localised level is at most occupied by one 
electron, and that the transition energy is appropriately 
rcnormalised due to the strong electron correlations. 

In the remainder of this Letter, the difference in the 
transport properties of the HIA and the loop correction 
in the non-magnetic case, is addressed. Therefore, let 
Tcr = r, and let ji = be the equilibrium chemical po- 
tential. In general, the renormalisation of the transition 
energies tends to increase the effective Coulomb repul- 
sion, e.g. U = A° 2S + (5A 2ff - A° - 5A a0 > U since 
A° 2S - A° = U and SA aQ < < 6A 2& . However, for 
U/r < 1 and U/T 3> 1, this does not cause any dra- 
matic deviation between the two approximations in the 
resulting transport properties. 

In contrast, by tuning the system into the regime 
1 < A° /r < A^/r (or A° /r < A^/r < -1) and 
0.5 < U/T < 1, for low temperatures, one finds that the 
current resulting from the HIA is weakly peaked for bias 
voltages such that either of the chemical potentials fiL/n. 
lies between the two transition energies, see inset of Fig. 
a) (dotted) . Consequently, there will be a region of 
a clear negative differential conductance (NDC) between 
the two conductance peaks corresponding to the tran- 
sition energies A a0 and A 2S , Fig. ^ a ) (dotted). This 
resonant peak is, however, not consistent with recent 
experimental data on single-level QDs, for example see 
Ref. [22l | . In the approximation with the loop correction 
(solid) , on the other hand, the differential conductance is 
positive in the whole bias voltage interval, in agreement 
with experiments. The current then shows a plateau, as 
expected when the one-particle state is resonant while 
the two-particle state is not. 

The difference of the two solutions can be understood 
as follows. The renormalisation of the transition energies 
yields A^o < A° and A 2S > A 29 . Hence, the state \a) 
begins to populate at lower bias voltages, Fig. ^b) (bold 



4 



0.06 

^0.04 
u 

■SZ- 

> 0.02 

;o 

■a n 



o 0.8 

\ 0.6 

\ 0.4 

7 - 0.2 




9 10 11 12 



— loop correction 
HIA 



(bold) - loop correction! ^) 



(faint) - HIA 




Support from Carl Trygger Foundation, Goran 
Gustafsson's Foundation, and Swedish Foundation of 
Strategic Research (SSF) is acknowledged. 



5 10 15 

bias voltage (mV) 



FIG. 1: (Colour online). Calculated transport characteristics 
of the single level QD. a) Differential conductance calculated 
in the HIA (dotted) and with the loop correction (solid) (the 
small ripples in the solid line is due to that dJ/dV is the nu- 
merical derivative of the calculated current). Inset: The cur- 
rent in the voltage interval between the conductance peaks, 
b) Occupation numbers No (dotted), N& (N-j = Ni) (solid), 
and N2 (dashed) in the HIA (faint) and with the loop correc- 
tion (bold). Here, {e , U, k B T}/Y = {1, 0.5, 0.02}. 



solid), than in the HIA, whereas the state \2) remains un- 
occupied for higher voltages (bold dashed). Therefore, 
the population in the one-particle state is able to sat- 
urate for bias voltages such that the two-particle state 
remains unoccupied. In the HIA, on the other hand, the 
one- and two-particle states compete about the avail- 
able population, since the transition energies lie closer 
to one another, which leads to an overpopulation of the 
one-particle state, Fig. \l\b) (faint solid). In general, a 
high population number in the one-particle state in com- 
bination with a significant reduction in the population 
number of the empty state |0), tends to reduce the tun- 
nelling probability of the one-particle state Hence, 
the overpopulation in |er) affects the transport properties 
of the system such that the current decreases. For larger 
bias voltages, however, the population in the one-particle 
state returns to its normal saturation value leading to an 
increasing current. 

In conclusion, an analytical formula for the current 
through a single level QD for arbitrary on-site Coulomb 
repulsion and couplings to the leads, was derived beyond 
standard mean-field theory (HIA) using a diagrammatic 
expansion of the QD GF in the strong coupling regime. 
The result removes the region of NDC, found in the HIA, 
in agreement with recent experimental data. 
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